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SYNOPSIS 


A number of coal cleaning processes ¥111011 are based 
on the float-sink principle are in use to clean a wide 
variety of coal changing both in quality and size consist. 

Por selecting a suitable process for a particular run— of- 
mine coal and forecasting the cleaning results, one needs 
performance criterion. Generally distribution (Tromp) curve 
which is obtained when the per cent of feed reporting to 
clean coal is plotted against density is used for this 
purpose. A formal model for coal cleaning processes that will 
result in distribution curve has been developed. The values 
of the parameters in the model are estimated using non-linear 
optimisation techniques namely Bo::- Comp lex method and 
Rosenhrock Hill Climb method. The variation of parameters 
with feed size gives a better insight into the performance 
of different coal cleaning equipments. 



CHAPTER 1 


IUTROSUCTI PIT 

Coal is the most abundant fossil fuel resource of 
energy. Raw coal is always associated with impurities. 

Gradual depletion of the high quality coal seams and 
increasing use of mechanised mining methods has led to the 
substantial increase in the impurity content of the raw 
coal, Frequently the raw coal is not acceptable for power- 
plants, metallurgical or other uses. A thorough cleaning 
of coal, to reduce the ash and sulphur impurities is now 
almost mandatory for coal utilisation. Because of the wide 
variations in the nature of coal, ash content, size fractions 
etc., coal cleaning is performed by various processes. For 
the selection of the optimum process and its monitoring, 
for a particular type and size of coal, some performance 
criterion is needed. 

Of the many performance criteria available for this 

purpose, distribution (Tromp) curve is most commonly used, 

\ 

In it the performance of various coal cleaning machines 
and the effect of control setting is given interms of error 
area, probable error etc. The mathematical form of the 
Tromp curve is identical to the cumulative distribution 
function. However, the interpretation of the two are 
different. In particular, the frequency function of Tromp 
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curve has no physical meaning. Many attempts have been made 
to describe Tromp curves interns of statistical cumulative 
distribution functions, including normal, log-normal and 
Weibull distributions. However, as of now a formal model 
for coal cleaning that will generate Tromp curve is not 
available. 

In this thesis, we present a mathematical model of 

the coal cleaning process that generates a separation (Trdmp) 

curve. The model is extended to incorporate the realistic 

case where the minimum densit3f of coal is greater than zero, 

maximum coal fraction floated is less than one and minimum 

coal fraction in sink is greater than zero. It is shown 

that the model’s separation curves are self preserving or 

self-similar interms of a dimensionless density, —3^ — , 

a 0.5 

inconformity with the results in literature. In the model, 
the sharpness of separation is quantified by a sharpness 
index n which is estimated using Box complex technique and 
Rosenbrock hill climb method. The variation of parameters 
with feed sizes in five different processes, dense medium 
cyclone, hydro cyclone, concentrating table, dense medium 
vessel and baum-jig has been determined and analyzed. 



CHAPTER 2 


LITERATURE REVIEW 
2*1 HEED EOR COAL CIEATOU; 

Coal is a mixture of combustible metamorphosed, plant 
remains that vary in both physical properties and chemical 
composition* A number of components aro present in the 
coal, depending on tho nature of the original plant materials 
and the degree of coalification. In addition tho coal 
coming from the pits and seams contain a significant amount 
of impurities, both ash forming mineral matter and sulphur 
contributing compounds. The origin, composition, properties 
and the distribution of impurities have boen discussed by 
many authors [ 1—3 ] . Shale, Kaolin, sulphide, carbonate 
and chloride group species are the major mineral impurities. 

A number of minor minerals, • quartz, feldspar, etc., are also 
present. The impurities are subdivided into 'inherent 
impurities' , which are structurally a part of the coal itself 
and cannot be separated from it by any mechanical/physical 
means, and 'segregated impurities', that exist as individual, 
discrete particles, and can bo removed by mechanical coal 
cleaning. 

Yancey and G-oer [ 4 ] havo summarized the ill effects 


of impurities: 
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1*' An economic loss is manifest in the handling, 
transportation, and storage of high-ash coal, for the inert 
impurities must be handled with the coal, yet contribute 
nothing in its use. 

2. In combustion, the ash-forming impurities not only 
dilute the combustible content of the coal and thus reduce 
its calorific value but also lower the effective capacity 
of the burning equipment, reduce the plant flexibility, and 
impose added expense in handling of ash. 

3. The fluctuations in ash content that characterize 
the high-ash coals render the control of combustion difficult. 

4. Excessive ash content impairs the efficiency of 
combustion and increases the erosion/ corrosion of the 
equipment . 

5. High ash content in the coal decreases the yield of 
coke interms of carbon content or thermal value, decreases 
the yield of by-products and gives a ■weaker coke and 
consequently higher percentage of breeze fraction. 

6. Use of high-ash coke in blast furnace leads to an 
increase in the consumption of coke, requires more limestone 
for flux, reduces furnace capacity, and makes the temperature 
control more difficult. 

7. Sulphur, the second impurity after ash is a dele- 
terious impurity in virtually all the metallurgical products, 
Eor this reason, high - sulphur coal is unsuited for smelting 
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or heat-treating processes. For example, in blast 
furnace operation, the penalities exerted by high-sulphur 
coke are added slag volume, increased coke consumption, 
and greater difficulty in furnace control. 

8. When coals containing more than about 0.05 par 
cent sodium chloride are coked, the refractory lining of 
the coke ovens suffers serious corrosion. 

The said ill effects can be minimized to a lesser or 
greater extend by coal cleaning which is defined as a 
physical and/or chemical process for removing high-ash 
impurities including pyritic sulphur from raw coal. The 
attributes of cleaned cool are greater uniformity in compo- 
sition and reduced impurity content. Zachar and Gilbert [5] 
have discussed the economics of cool preparation in detail. 
According to them, the removal of free impurities from raw- 
coal to upgrade the product inorder to make it saleable is 
an expensive process employed only through necessity: 

1. to increase its net income per ton of product, and 

2. to provide a steady outlet for its products. 

Mechanization and the use of continuous mining 

equipment invariably produce raw coal containing more im- 
purities. Coal preparation therefore permits rapid, full- 
seam mining. Also the amount of impurities in fine coal 
tends to be high, which again oalls for appropriate coal 
cleaning facilities. 
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2,2 COAL CLEANING PROCESSES 

Most of the coal cleaning processes depend on the 
density difference between coal and its associated impurities. 
Coal ranges in apparent density from 1,2 to 1*7, depending on 
its rank, moisture content, and percentage of ash. Shale, 
clay, and sandstone range from about 2,0 to 2,6 in density, 
depending on their degree of purity; pyrite ranges from 2,4 
to 4,9; and cal cite and gypsum have a, density of 2,7 and 2,5, 
respectively. All removable impurities associated with 
coals are higher in density than the co.a.1, The two major 
difficulties, that have influenced the development of a 
number of coal cleaning processes are: 

1, The raw coal contains middlings of a wide range 
of intermediate ash contents; and 

2, The relative behaviour of the coal particles of 
different sizes. 

The classification of the coal cleaning equipments 
and the principles underlying them have been dealt by many 
authors [6-9]. A simple classification of the more important 
wet, gravity concentration methods/ equipments is: 

1, Dense medium vessels 

2, Jig washers 

5. Dense medium cyclones 

4. Hydro cyclone 

5, Concentrating Tables, 
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Dense medium vessels and jig washers arc normally 
used for the concentration of coarse coal. Palowitch and 
Dourbrouck [9a] have described the dense medium process, 
the separator design, the various kinds of equipment available 
and 'the advantages of dense medium processes over other coal 
cleaning processes. A major portion of the cools is washed 
by dense medium processes, because of the following chara- 
cteristics of this process: 

1. Ability to handle a wide range of sizes 

2, Ability to make sharp separations at any density 
within the normal range, even in the presence of 
high percentages of the feed in the range of +0.1 
density units. 

A major disadvantage is tho,t the fine sizes cannot be 
treated, because they decrease the density of the med^tm and 
cannot be recovered from the sand (suspension material) by 
screening. 

The unit operation of jigging, the different types of 
jigs used and the design of Baum-jig have been given by 
Lovell [9b]. Jigs are employed to clean either sized or 
Unsized coal; in a few instances they have boon used to 
treat in one operation all sizes from dust upto 6 or even 8 
inches. When such a wide size range is employed, the cleaned 
product may meet the market requirements, but the efficiency 
with which some of the individual sizes are treated, is 
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necessarily limited. Jig cleaning efficiency tends to 
decrease -with, the smaller sizes and tends to clean best at 
higher densities as size decreases. When a very broad size 
range is cleaned by a jig, the efficiency of separation 
decreases 'with a decrease in average size. According to 
Yancey and G-eer, the recovery efficiency was nearly as good 
as attained with heavy media jirocess. 

The cyclones and concentrating tables are mainly used 
for cleaning of fine coal. Sokaski, G-eer and Yancey [9c] 
have described the theory, flowsheet, description and per- 
formance' of dense medium cyclones and hydrocyclone. In general, 
the dense medium C 3 '-clone effects a sharper separation between 
coal and impurity than can bo obtained in other typos pf 
cleaners handling the same size range. 

As of now, cyclones are used to treat a comparatively 
narrow size range, typically 3/8 or 1/4 inch to 1/2 mm. 
Therefore, little information on how various sizes behave 
in the cyclone is available. It would bo however expected 
that coarse coal would bo separated more sharply than fine 
coal. When cleaning is performed at a high density setting, 
the superior sharpness of separation accomplished in the 
cyclone does not greatly improve the recovery efficiency and 
the economics of the operation. However, if the separation 
must be carried out at a low density, cyclone alone can 
produce washed coal of required sulphur or ash content, at 
higher efficiency. 
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The separations obtained in a hydro cyclone are not 
nearly as sharp as those that are characteristic of the 
dense medium cyclone. Consequently hydro cyclones are not 
applicable for difficiilt to clean coals or separation at 
low density. The hydrocyclonc could especially be Suited for 
treating minus 28-mosh coal, more so if the coal is not 
amenable to flotation. The hydrocyclonc is reported to bo 
superior to flotation for lowering the sulphur content of the 
washed coal, if iron pyrito is present in the food. 

Deurbrouck and Palowitch [9d] have described the 
principle of tabling processes, opera/ting conditions and 
their performance characteristics. The concentrating table 
cleans the fine particles much more efficiently than either 
hydro cyclones or the fine coal jig (feldspar jig). However, 
minus 200-mesh material was found not responding to upgrading 
by tabling. 

Using the flow-sheet of Bhojudih washery, located in 
the eo„stcrn part of Jharia coal field, G-okale and Rao [10] 
have illustrated, how the different cleaning processes are 
combined in a coal washery. Here coal is crushed to minus 
3-inch size and screened on 1-inch screen. The minus 3-inch 
plus 1— inch fraction is cleaned by heavy media separation and 
the minus 1-inch fraction is washed in baum jigs. The sinks 
of heavy media and the jigs ore then crushed to minus 1/2 — 
inch and washed in dense media cyclones. 
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Based on the data for various schemes of washing in 
vogue at different washerics as compiled Toy Ghakravart i 
and Moritra, G-okale and Rao observe the following points: 

1, The accuracy of separation at any predetermined 
density within the range of 1.25 - 2.50 and continuous 
maintenance of preselected value of density within + 0.005 
of its value, are the salient features of dense medium 
processes. The additional qualifications of heavy media 
separation arc its ability to separate at one or more pre- 
determined densities and its capacity to treat a larger 
size range of coal, besides its low cost, flexibility, and 
accuracy of control. 

2. Baum-jigs arc employed where initial separation of 
shales from coals is required* The separation is effected 
at a density between 1,70 and 1.80. 

5. High ash content and inferior coals require the 
flotation treatment, 

4. The dense medium cyclones can take care of the finer 

size (less than 1.25 ems in size) coals of difficult washing 

<**• 

characteristics, with an optimum recovery* 

2.3 THE PEREOEHAHCB CRITERIA 

It is evident that the selection of the best process 
flow for a given run-of— mine coal requires a careful study 
wherein many conflicting factors must be weighed. The cost 
of a detailed investigation is well repaid capital investment, 
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in higher recoveries, In flexibility, and in the ease of 
operation and maintenance. In order to compare the different 
cleaning operations, the performance of coal cleaning equip- 
ment or kinds of coal, and to forecast the cleaning results 
we need a suitable performance criterion. 

2.4 TYPES 0? PEREQBMAECB CRITERIA 

Geer and Yancey fll] have classified the various 
performance criteria according to the degree to which they are 
dependent on the density composition of the raw coal, namely 
directly dependent, indirectly dependent and independent 
criteria. 

The yield and ash content of the washed coal are two 
performance factors of direct and immediate interest. These 
two factors are directly dependent on thedensity composition 
of the raw coal treated and the functioning of the equipment. 
The continuous variation in the amount of impurity in the 
raw coal makes it impossible to tell whether a reduction in 
yield or an increase in the ash content results from any 
malfunction or mi sad justment of the cleaning equipment. As 
such these performance factors are not suitable for comparing 
the performance of different plants, because no two plant 
feeds are the same. 

The important performance criteria, that are indirectly 
dependent on the density composition of raw coal are; 
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1. Efficiency formulas 

2. Misplaced material concept, and 

3. Ash error and Yield error. 

The well known efficiency formula given by Eraser 
and Yancey is. 


Efficiency = 


Yield of washed coal 

Yield of float coal of same ash 
content 


x 100 


This type of efficiency formula has a drawback of 
giving a value of 100, when no impurity is removed. In 
pn-eumatic cleaning equipments impurities finer than 28-mesh 
size are not removed. So fo.r such fine size raw coal, the 
efficiency becomes 100. To check such drawbacks, ash reduction 
must be considered along with efficiency. Satisfactory 
performance then necessiates both high efficiency and the 
required ash content in the washed coal. Also to account 
for the coal that is freed by degradation during the washing 
process, one has to use calculated feed in the efficiency 
formulas. 

'Misplaced material' is the sum of the sink in the 
washed coal and the float in the refuse, generally expressed 
as a percentage of the raw coal. Although the displaced 
material* concept is simple to understand, it does not take 
into account the quality of the misplaced material. It does 
not relate directly to the all important factors of yield and 


ash content 
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Ash. error is the numerical difference "between the 

/ 

actual and theoretical ash contents of washed coal at the 
yield of washed coal obtained. It is a direct measurement 
of tine impairment on account of ash content. Yield error 
is the difference between the yield of coal actually obtained 
and the theoretical yield at the ash content of the crashed 
coal. It expresses the loss in j^ield. The authors [ll] 
have given the general relationship between efficiency, mis- 
placed material, ash error and yield error, in the form of 
a tabulated data. The restrictions on the applicability of 
these directly and indirectly dependent criteria are stated as 
follows: 

’Equipments treating coals of substantially similar 
density compositions or operating at substantially same densities 
of separation alone can be compared on the basis of these 
crit eria . * 

2.5 PIS TRI BUT I Oh OR SEPARATION ( TRQHP ) CURVE 

TJie type of separation between coal and impurity 
obtained in actual practice is illustrated graphically in 
Eig. 2.1* Coal of low density and impurity of high density 
report largely or entirely to their proper products. As the 
densitjr of separation is approached, more and more material 
reports to the improper product. The infinitesimal increment 
of density that is divided equally between the washed coal 
and refuse is defined as the density of separation. 
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Tromp observed that the curve in Fig* 2,1, which has 
the shape of a crude Gaussian error distribution, is 
governed by the laws of probability. Also, if the percentage 
of recoveries of washed coal are plotted against the mean 
densities of the density fraction, the type of curve in 
Figure 2,2 is obtained. The curve in Figure 2,2 has many 
names such as Tromp curve, distribution curve, separation 
and recovery curve. The distribution curve in Figure 2.2 
can also be obtained by rotating the left-hand limb of curve 
in Figure 2.1. by 180°, about the 50 percentage ordinate. 

When the distribution curve is plotted on either log 
probability or arithmatic probability grid, one will occa- 
sionally get a straight line plot. The advantages inherent 
in a straight-line plot of this kind are: 

1. The slope of the line is a measure of sharpness of 
separation. Steeper the slope, sharper is the classification. 

2, The slope and the specific gravity of separation 
combine to characterize the complete curve. 

However, in practice, many distribution curves do not 
exhibit a straight line throughout the full range on log 
normal or normal plot. 

As stressed by Gottfried [12] , distribution curve is a 
plot of the weight per cent of feed reporting to clean coal, 
f(d), as a function of density, d. The yield of clean coal, 

Y is, 
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f 1 

) f(d) m(d) dd 



J m(d) dd 

6 


Itfhere m(d) represents the feed rate of material having 
density d and d-^ represents the highest density value of the 
feed material* 

The density of the material in the raw feed that is 

divided equally "between the clean coal and refuse is defined 

as the density of separation, d . In otherwords, it is the 

s 

value of density corresponding to f(d) = 50. By making appro- 
priate physical adjustments in a coal washing machine, the 
value of the density of separation d s can "be increased or 
decreased, thus shifting the entire distribution curve to the 
right or left. 

Gottfried and Jacobsen [l3] have discussed the distri- 
bution curve for a theoretically perfect separation and an 
actual separation case as in Figure 2.2. The shaded area, 
called error area,, in. Figure 2.2 is a measure of the difference 
between the actual and theoretically perfect separation. The 
error area approaches zero as the actual distribution curve 
approaches the theoretical, perfect classification. Error 
area is dimensionless and ranges from as low on about 10 for 
good dense medium operation to over 100 for some inferior 
table and jig operation. 
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Another measure of the distribution curve called 
‘Probable error’ also represent the sharpness with which the 
coal and impurity are separated. The quantity is equal to 
one-half the density difference between the 25 and 75 per cent 
ordinate values on the curve. The steeper the distribution 
curve, tho lower the probable error. Pigure 2,2 illustrates 
the relationship between probable error and the slope of the 
distribution curve. 

Geer and Yancey [ll] have compared the distribution 
data for four dense medium plants and concluded that within 
reasonable limits, the separation between coal and impurity is 
independent of the density composition of the feed. For other 
processes also the independent nature of the distribution curve 
has been established. 

The size composition of the raw coal influences the 
distribution curve. In dense medium vessel and dense medium 
cyclone, coarse sizes are separated somewhat more sharply than 
tho fine sizes. In jigs and tables this dependence is even 
more pronounced. The distribution curves for the coarse, 
intermediate and fine sizes of coed, cleaned on a concentrating 
table have been given by Geer and Yancey [ll] • 

2.6 GENERALIZED (SELF PRESERVING) DISTRIBUTION CURVE: 

Por a given coal-cleaning equipment and a given size 
distribution of coal, different distribution curves as in 
Figure 2.3 are obtained for different settings of density of 
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separation* The sliape of these different distribution 
curves are similar, because of the following factors. 

For a given size distribution of coal and machine 
1* Tromp curve is independent of density distribution of 
coal and 

2 . Tromp curve is independent of density of separation, ^ 
setting. 

A method for combining those different distribution 
curve into one generalized or self preserving curve, which 
will apply for any density of separation within the range of 
original data, is presented by G-ottfried and Jacobsen [l3]« 
Here tho weight per cent of feed reporting to clean coal; f(d), 
is plotted against the reduced density, d. The reduced dimen- 
sionless density of separation is defined as the ratio of 
density to the density of separation. 


l. e. 


-d = 


d 



The resulting distribution curve is independent of the 
density of separation. They have shown, how the several 
conventional distribution curves of given coal-cleaning equip- 
ment treating particular size distribution of coal, can be 
represented by a single generalized distribution curve gra- 
phically. The significance of tho generalized distribution 
curve is that it can bo used with any density of separation, 
not only the density of separation corresponding to tho given 
data. 
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The effect of size on the generalized distribution 
curve, the uses and range of applicability of generalized 
distribution curve have been well dealt by Gottfried and 
Jacobsen [l 3]« Figure 2.3 shows conventional distribution 
curves and Figure 2*4, the resulting generalized distribution 
curve , 

2.7 MATHEMATICAL FORK OF DISTRIBUTION CTJFVE 

The various attenpts at mathematical representation 
of distribution curve have been dealt by Gottfried [12]. 

Tromp recognized the similarity between the distribution curve 
obtained for a coal cleaning device and a normal (Gaussian) 
distribution function. This normal distribution function is 
not suitable because. 

1, The cumulative normal distribution function cannot 
be expressed in a closed form (that is, in terms of a simple 
algebric formula). This is also true for log normal distri- 
bution. 

2. Many coal -cleaning devices exhibit asymmetric 
distribution curves, whereas the normal distribution is 
perfectly symmetrical. 

Gottfried [12 ] has suggested that the weibull function, 
an exponential type distribution function may turn out to be 
a better candidate for a mathematical representation of the 
generalized distribution curves. 

1. Weibull function has several interesting character- 
istics, including a closed-form representation for the 
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cumulative (integrated) function. 

2. With, proper choice of parameters, it can provide 
a reasonably near closed-form approximation to the cumulative 
normal distribution function. 

3 » It can also be used to represent asymmetric 
distribution curves. 

'The cumulative (integrated) distribution form of 
Weibull frequency function is, 

f (d) = 100 exp [-(d - d 1 ) a /b] (2.1) 

Where d is the reduced density, f(d) is the per cent 
of feed of reduced density d reporting to clean coal, a and b 
are constants. Equation (2.1) will apprpximate the genera- 
lized distribution curve. Since in actual practice of coal- 
cleaning, the maximum floated coal may not be 100 per cent 
but something like 98 or 99 per cent and the minimum floated 
may not be zero but as high as 30 or 40 per cent, equation 
(2.1) is modified as follows [l2] * 

f(d) = 100 [f 0 +c exp -(d -d 1 ) a /b ] (2*2) 

Where f Q +c is always less than or equal to 1. Applying 
the condition that f(d) = 50, when d = 1, 

c l/a 

d-, = 1 - [b log ] ( 2.3) 

1 e 0.5- f Q ~ 

Hence, the number of independent parameters is '4* 

that is f , a, b and c. 
o 7 1 
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Equation (2.2) can Toe used to represent the generalized 
distribution curve, provided a,b,c and f are determined- for 
a given set of distribution data. By nonlinear regression 
technique Gottfried [12] computed the values of a,b,c and f 
and tabulated them for several different cleaning vessels 
and a number of feed size fractions. 

Although the fitted weibull functions represent the 
experimental data reasonably well, it has a number of dis- 
advantages, 

1. Many of the curves broke too sharply at the extreme 
points, that is at the knees of the curves, 

2. The fitted curve often fall short of spanning the 
complete range namely 0-100 per cent. The disparity between 
the analytical and the tabular forms of the distribution data 
is quite severe in some curves. 

3. The choice of the Weibull distribution is based on 
an adhoc assumption, and it is difficult to relate its 4 para 
meters to the underlying physical process of coal cleaning. 
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CHATTER 5 


OBJECTIVE 

Prom the previous chapter on the literature review 
of this subject, we conclude that there is at present no 
formal model on coal cleaning processes that results in 
Tromp curve. Our objective in this work is to develop a 
model for a class of coal cleaning processes from a set of 
phenomenological concepts of the cleaning machines. We show 
that the model generates a distribution curve which is self- 
similar. That is, the separation is independent of the 
density composition of the raw coal feed and reduced density 
of separation, within reasonable limits. Further, the model 
is tested on computer using the rather extensive data provided 
by Geer and Yancey [ll] and Gottfried and Jacobsen [l3 3 • The 
sharpness index is estimated and analysed for five kinds of 
coal washing machines and a broad range of coal size fractions. 



CHAP P Ml 4 


MOP EL 

Most coal cleaning oqxd.pnon.ts employ a float-sink 
principle, based on the density difference between coal and 
associated impurities. Thus a typical piece of coal cleaning 
equipment will separate a given feed into a low density clean 
coal product and a high density refuse. The schematic diagram 
of a float— sink process is given in Eigurc 4«2-« 


m " (a) fio=.t 



sink 


Ei,?ure 4.1 

let m* (d) be the mass flow rate of the raw coa.1 of 
density d that is fed into the cleaning equipment, m* 1 (d) be 
the mass flow rate of the washed coal, m ,n (d) be the mass 
flow rate of the impurity and m(d) be the steady state mass 
of coal in the suspension in the machine* 

Since the rate of flotation depends on (1) mass of 
coal that is present in the suspension (2) density of coal in 
the suspension and machine characteristics, we can incite. 
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Rato of float a m(d) 

.. a 

a d 

Therefore, Rato of float = K* d a m(d) 

Where a is some exponent and E*' is the machine 
constant. 

Similar rate of sinking of impure material is, 


or. 


i« e, 
i.e. 


Rate of sink 


Rate of sink 


a 

a 


m(d) 


= E" m(d) d 


Where E* is also a machine constant. 

At steady state, by mass balance, 

Rate of feed = Rate of float + Rate of sink 


m* (d) = IV d a m(d) + E 1 ' d u m(d) 
m’(d) - E* d a m(d) - E* d d m(d) 


0 


Therefore m(d) = 


m* (d) 


feed 


E‘d A + E” d b 


Since m M ( < 3-)-f2.oet = ^ 

Substitution of equation (4.2) gives, 

n ’ ^ d ^feed 

m'* (d)^ . = v 

float K , d - + E ji 


E* d‘ 


a 


(4.1) 

(4.2) 

(4.3) 

(4.4) 


Hass fraction recovered is, 

m " Afloat 
B '^ d) fooa 


f(a) = 


(4.5) 
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Substituting the values of m 11 (d)-, 

xloci 


f (a) - . __i_ 

E' d a + K" d D 1 + E d n 

where, 

K * 1 

E = and n = b~a 

E»' 


(4.6) 

(4.7) 


\Ic know at the split point, that is at the density 
of separation d^ 50 per cent or 0,5 fraction will be 
recovered. That is f(d Q = 0,5 

Substituting this value of f(d Q in the equation 
(4.6) wo get, 


0.5 = 


Therefore E = 


1 + E d 


n 


0.5 


,n 


(4.8) 

(4.9) 


0.5 


Substituting this value of E in equation (4.6) leads to 


f(d) = 


1 +(• 


_d_ 

^ 0.5 


) n 


( 4 . 10 ) 


let a dimensionless density d be defined as d = v — , 

a 0.5 

then equation (4.10) becomes, 


f (3) = (4.H) 

1 + d 1 

This is the basic model of floc/t-sink processes. Here 
f(d) is the fraction that is recovered, which is a function 
of some dimensionless densitjr d and n is the sharpness index. 

If the density of coal considered is zero, the reduced 
density d is zero and the fraction recovered becomes Ij that 
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Is all the f e„e.d will "be floated. If the density considered 
is infinity, d is also infinity and the fraction recovered 
becomes zeroj that is nothing of infinite density mil be 
floa.ted. If the density of coal under consideration (d) is 

.the same as that of density of separation, that is then 

— 5 

d = -5 — ^ = 1, and the fraction recovered will be 0.5. 

a 0.5 

Using the basic equation (4.11) of float-sink process, 
the nature of the distribution curves for various values of 
the sharpness index, n are determined, using tho programme in 
Appendix 3. The computed values of f(d) for different sharp- 
ness index, n are given in the same appendix. The same results 
have also been shorn graphically in Figure 4.2. As the sharp- 
ness index, n, value increases, the distribution curve becomes 
steeper and approaches tho ideal or theoretically perfect 
separation. 

Bote that n, the index of separation, is always greater 


than zero. For 

the case, 

n is infinity; 

... 

(i) d < 

d 0.5 ' 

d 4 1 and 

f(d) = 1 

(ii) d > 

d 0.5 ’ 

d 1 and 

f(d) = 0 

(iii) d = 

d 0.5 ' 

d = 1 and 

f(d) = 0.5 


This is the ideal case for perfect separation shorn 
in Figure 2,2. 

The above basic model is based on; 

(1) Minimum density is zero, 

(2) Maximum float is l f and 
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(3) Minimum float is zero. 


But in actual practice the minimum density may have 
some value, d 1( greater than zero, the maximum float f , 
may he less than 1, and the minimum float f^, greater than 
zero (See Figure 4.4). 

Wo now proceed to modify our model to extend its 
validity. First, still assuming f^ = 1, and f = 0, hut 
d-^ > 0, we developed a 2-pa.rametcr model in the following 
manner. 


We modify equation (4.1) as, 

m‘( d )fe ed “ K’ (d- d i) a m(d)-K’« m(d) (d-d-^ = 0 

( 4 . 12 ) 

That is instead of considering d alone, we consider (d~d-^) 
Hence equation (4.6) becomes 


f(cL) = 


(i+£(a - a 1 r) 

When f(d) = 0,5 and d - d Q ^ 
0.5 = 


(l+KCd^j-d^ 11 ) 

and equation (4*10) becomes 


(4.13) 


(4.14) 


f(d) 


f(3) 


l + ( 


d-d. 




d 0.5- d l J 



(4.15) 


or 


( 4 . 16 ) 
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Let us. note, in this 2-paraneter model also_fraction 
recovered f(d) is a function of reduced- density, cl and not 
of density d. Here also f(d) = 1 when d = d^ and f(l) = 0.5 
and f( ao ) = 0. 

Based on the 2-parameter model equation (4.16), the 
effect of variation of d^, (the minimum reduced density of 
coal that will result in 100 per cent float) for a specific 
sharpness index, n is studied. The computed results for the 
above effect have been given in Appendix 5. The same results 
have been shown graphically in- -Figure 4*5. Prom Pigure 4.3 
it is evident that higher the d^ value, higher is the 
separation or cleaning of coal. 

How we come to the general case of a 4-parameter model, 
talcing into account f < 1, f^ > 0 and d^ > 0. Por this 
purpose we make a transformation as follows: 



Here also f = 1 when f= f u and T = 0 when 
Equation (4.15) becomes 


(4.17) 


f= 


f - f. 


f u“ f l l+Etd-d^ 11 
At the split point. 


0.5..,- f 


1 = 


f u - f i 1+E(a 0>5 - a 1 ) n 


(4.18) 


(4.19) 


* - r ( f n ) 1" 1] ^ 

K - u 0.5 - (d 

1 - ^ a 0.5 


d l> 


n 


or 


(4.20) 
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Therefore equation (4.13) becomes 


f - f 1 

f - f 

x u X 1 


f ~0.5 r 

< s 


1 - 


4i) 

d 


( 4 . 21 ) 


n 


or 


f = f(d) 


f l + 


f --fi 
u 1 


f - 0.5 v 

i+(A^ ~) (^ 1 ) n 

l-d. 


‘0.5- 


( 4 . 22 ) 


Here also we note, 
f(di) = f u 
f('-xi) = f x 
f(l ) = 0.5 

and f(d) depends on d and not on d, that is the fraction 
recovered is a function of reduced density is never violated 
in the setting up of 2- and 4-parameter models. 

A representation of 4 parameter model is given in 


Figure 4.4 
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CHAPTER 5 


methodology 

In the Chapter 4, we have developed a 2— parameter 
and a 4-parameter model for the analysis of coal cleaning 
operation. The 2-parameter model is 


f(d) = 


-1 


1 +( - - ^ ~l) n 


1-d. 


( 4 . 16 ) 


where d-^ and n are the two parameters that must be deter min ed 
from the experimental data. The 4-parameter model is. 


f(d) = t ± + 



1 +( 


ye- 5 

0.5-f 1 



n 


( 4 . 22 ) 


where d^, n, f^ and f^ are the four parameters to be determined. 

The parameters in our models were estimated from the 
rather, extensive data provided by Gottfried and Jacobsen [l3], 
where d and f(d) values for various sizes ranges of coal 
processed in different plants are given. The breakdown of the 
data are given in the following Table 5«1. 

Our objective is to estimate the values of 2 and 4 
parameters in the two models from this data that will represent 
the physical situation, such that f(d) calculated using 2- 
or 4-parametor model should be as * close* to f(d) experimental 
(which is available in data) as possible. 
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TABL E 5.1s GEHE HAIIZED D ISTRIBHT10II DATA [l3. 


S .Eo 

1 . 

2 . 


3. 


4. 



Equipment 


Peed sizes, in. 


Dense medium vessel 


Baum jig 


Dense medium 
cyclone 


Hydro cyclone 


Concentrating 

table 


6x4, 4x2, 2x1, 1x1/ 2, l/2xl/4 
and composite feed (6xl/4). 

6x3,3xl| , l|xl/2, 1/2x1/ 4, 

1/4x8 mesh, 8x14 mesh, 14x68 mesh 
and composite feed (6 x 48 mesh). 

3/ 4x1/2, 1/ 2x3/8, 3/8xl/4, l/4x8mesh, 
8x14 mesh, 14x28 mesh and composite 
feed (3/4 x 28 mesh). 

l/4x4mesh, 4x8mesh, 8x4mesh, 

14x28 mesh, 28x48 mesh, 48x100 
mesh, 100x200 mesh and composite 
feec ( 1/ 4 x 200 mesh). 

3/ 8x1/4, 1/4x8 mesh, 8x14 mesh, 

14x28 mesh, 28x48 mesh, 48x100 mesh, 
100x200 mesh and composite feed 
(3/8 x 200 mesh). 



38 


Following the standard proccduro For parameter esti- 
mation, the objective function can be written as 
n 

(f^(d) - f ± (d )f (5.1) 

i=l 

where i = 1,2,... n, and n is the number of experimental points 
considered. By minimizing the above objective function, we 
can produce the best estimate of 2 and 4 parameters in our 
models. The optimization techniques used for the above 
purpose wore 

1. Box complex technique [l4,15] 

2. Rosenbrock hill climb method [16,17] 

The two methods are discussed, with flow chart and 
computer program, in Appendix 1 and Appendix 2 respectively. 

5.1 Selection of 4-Parameter Model 

For a particular coal cleaning vessel and a given size 
range of coal, both 2-parameter and 4-parameter models were 
tried. Almost all the trials indicated that the 4-parameter 
model reproduces the experimental results better than ohe 2 — 
parameter model. Example of one such case is given in Table 
(5,2). So we selected the 4-parameter model for future 
analysis in this study. 

5.2 Selection of Best Weighted Objective,. Fun ction 

The f(d) data which ranges from 0.0 per cent to 100.0 
per cent may not be very accurate at both extreme ends. To 




coupons ci *oG for fiiis, fhc weighted objective functions were 
also tried, nanoly 
n 

$ w (f ± (d) - f ± (d)]p (5.2) 

i=l 


where the ¥ is the weight given as 

¥ = (50 - Abs (f ± (d) - 50))/50 

This weight will give a value of 1 when 11(d) = 50, 
and a value of zero at f£d) = 0 and 100. That is, more 
weightage is given to the central portion of the data than 
the end points in the two tails. In some trials the extreme 
points were heavily weighted in order to force fit the curve 
in the tail region. 

The weightage for the extreme points can be given by 
just replacing ¥ by (l-¥) in the objective function. For 
uniform weightage throughout the experimental data points, 
value of unity is given for ¥. 

The extent to which the model reproduces the experimental 
data was determined interns of absolute error per data point 


and normalized error per data point: 


(f ± (d) - f^d)) 2 


Absolute error = 


Normalized error 


^ (f ± (d) - f-^d)) 2 
L=1 t ± (d) • 
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Many trials conducted in this manner showed that the 
objective function weighted uniformly gives the best re- 
production of experimental results. So, we selected uniform 
weight, that is no weight is imposed on any part of the 
distribution curve. Example of uniform weight, weighted at 
mid portion and weighted at extreme points were given in 
Table 5.3. 

It was also determined that the absolute error criterion 
in general gave better results than normalized error criterion. 
Example to the effect is given in Table 5.4. 

Ror every data set, estimation of parameter was carried 
out for 0-100 per cent, 1-99 per cent and 2-98 per cent range 
of f(d). As we expected, the 2-98 per cent range gave best 
reproduction. Example of this is given in Table 5.5. 

5.3 SoleotioKL-Of. Roscnbrock Mill Climb Method 

Eor the 4-parametcr model, and uniformly weighted 
objective function, both Box Complex technique, and Rosenbrock 
Hill Climbing I-lothod were tried for each data set. It was 
found that in some cases both methods give identical results. 

In other cases however, for box-complex method to be success- 
ful, (i) the reflection coefficient, the expansion co- 
efficient and the contraction co-efficient had to be properly 
selected, (ii) moreover the success of the method was quite 
sensitive to the choice of the initial guess points, and 
(iii) also box— complex technique took more computer time than 
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tlic Rosenbrock bill climb procedure. Consequently, Rosenbrock 
bill cliiab procedure was used for estimating Hie 4— parameters 
in our model* An example to tbis effect is given in Table 5*6, 

The maximum and minimum values of the 4-parameters 
should be provided for the estimation of the 4-parameters, 

Prom the studies on nature of distribution curves for various 
sharpness index, n as in Appendix n around 50 is found to 
give experimental conditions. So the range given for n is 
0-100. Since d^ can be any value between 0 to minimum 
reduced density experimentally found, the range given for 
d-^ is 0- d^ when d^ is the minimum reduced density in data 
for which the 4 parameters should be estimated. Also to give 
a wider choir e, f u (d) is given a range between 50 - 100 per 
cent o-nd fq(l) between 0 - 50 per cent. 

Since the optimization technique used in the estimation 
of 4 parameters is heuristic sometimes ire got too high n 
values and too low d^_ values. Such unrealistic values are 
discarded and reasonable values of the 4 parameters were 
selected and were listed out in Appendix 4. 
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CHAPTER 6 


RESULTS 

The theoretical model developed in Chapter 4 for 
coal cleaning processes was tested on computer using the 
rather extensive data provided by Gottfried and Jacobsen [13]. 
Based on the trials reported in Chapter 5 on methodology, it 
was decided to use Rosenbrock hill climb method for the esti- 
mation of 4 parameters in the model. Erom the estimated 
values of the parameters for a given size of raw coal feed 
in a coal cleaning equipment, the absolute error, standard 
deviation and the range of washed coal percentage, f(d) 
covered were computed.' The mathematical representation of 
the distribution curves suggested recently by Gottfried [12], 
by Weibull function were also used to find similar absolute 
error, standard deviation and the range of f(d) covered for 
the same sets of experimental data. The comparative results 
of both Gottfried equation and our model for all five types 
of coal cleaning equipment and various sizes of coal feed are 
given in Appendix 5» Figures 6.1 — 6*5 show the comparison 
between the two models. 

It is evident from the above results that our model 
covers a wider range of washed coal percentage f(d) than it 
is possible by Gottfried equation. Also the absolute error 
and standard deviation are far less in our model. Invariably 
even within a given range of applicability our model is found 
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to be better in reproducing the distribution curves. It is 
to bo noted that in our model, unlike Gottfried equation, 
the distribution curves do not break sharply at the extreme 
points, and the disparity between the predicted end experi- 
mental points in the distributions is never very large in 
our model. ■ Moreover in our model the 4 parameters have 
physical significance, namely, 

f u (d) - Maximum washed coal percentage that can be 
expected for a particular size range of raw coal feed in a 
given coal cleaning equipment, 

f-^(d) - the minimum coal percentage that will be in 
the sink 

-•n - the sharpness of separation 
and d-j^ - the minimu m reduced density of the coal that will 
completely float in the process. 

The four parameters of our coal cleaning process 
model have been estimated for all five types of coal cleaning 
equipment and all sizes of coal in the given da,ta. The values 
of these parameters, the absolute error and standard deviation 
are given in Appendix 4» Figures 6.6 — 6.11 show the effect 
of feed size on sharpness index, in different coal cleaning 
equipments. 
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CHAPTER 7 


ANALYS IS 

7.1 Dense Medium Vessel 

The sharpness index n varies from 12.3 to 13.3 for 

feeds in size ranges 1/2” x 1/4” to 4” x2* as given in 

Figure 6.6. Raw coal feed of sizes lower than 1/4* 1 and 

* * 

greater than 4** have lower sharpness index. The composite 
feed of size 1/4” x 6” has a rather high sharpness in dex 

* * j 4 

showing the ability of dense medium vessel to handle wide 
size range of raw coal feed. On the other hand, the inability 
of dense medium vessels to handle fine size is evident from 
the droping nature of sharpness index curve below 1/2* size. 
Those conclusions regarding dense medium vessels are in 
conformity with the literature on dense medium vessel per- 
formance [9a], 

f^(d) - the maximum washed coal percentage that can be 
expected varies from 96.5 to S9.9 and f-^(d) - the minimum coal 
in sink varies from 0.11 x 10“^ to 3*03 per cent. The value 
of d^ , the limiting reduced density varies from 0,75 to 0.94 
.indicating tha/t there is no difficulty in cleaning, 

7.2 Baum Jig 

In agreement with the literature [9b], the sharpness 
index Increases with increase in feed size as given in 
Figure 6.7. It varies from 2.8 to 20.6 over a Size range 
14 x 28 mesh to 6* x3 M and for a compqsite feed of size 
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\ varied from °* 76 to 0«9, even difficult to wash coals 
can he handled by this process, 

7.4 Hydro cyclone 

The variation of sharpness index n with feed size 
is given in Figure 6.10. Here the sharpness index varies 
from 6,5 to 10,22, and is lower than for dense medium cyclone. 
The speciality of hydro cyclone in treating relatively fine 
14 x 28 mesh size feed is evident by the parameter value of 
n = 10.2. Raw coal of feed size greater or less than this 
28 mesh size has lower sharpness index. The composite feed 
of size 1/4* x 200 mesh has a very low value i.e. n = 2.39. 

•'ll 

Clearly hydrocyclone is found to he unsuitable for broad 
size range of raw coal feed. 

Here f (d) varied from 92.6 to 99.88 indicating some- 
what less recovery than dense medium cyclone. The f^(d) 
varied from 0.005 to 10.9 indicating that under certain 
conditions more than 10 per cent of valuable coal may be 
lost to sink by this process. The variation of d^ from 0.11 
to '0.5 show that hydrocyclones are not suitable for difficult 
to clean coals. 

7.5 Concentrating Table 

As shown in Figure 6.11, the sharpness index, n 
fluctuates sharply from 4.56 to 10.92. The underlying 
principles of concentrating table are somewhat different from 
the equipment considered earlier. It is not clear if our 
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model has a meaningful relevance to concentrating "tables* 

In any case, here f^(d) varies from 97,5 to 99,99 
indicating higher recoveries of wasted coal, however, f^d) 
values sometimes rise as high as 18, 4 , showing the m— 
reliable nature of this process, 

7,6 G-en eralized Probable Error : 

The conventional probable error is given in Pig, 2, 2, 

In the case of generalized distribution curve as in Figure 2,4, 

we will have generalized probable error (GE ) which can be 

P 

simply written as 

6 *, = 0.5 [a 25 - J ?5 ] 

where d^ is the reduced density at which 25 per cent 
of feed coal will report to clean coal product and d^- is the 
reduced density at which 75 per cent of feed coal will report 
to clean coal product. 

Using the numerical values of the 4 parameters of 
our model, generalised probable error for all feed size • 
fractions were calculated and tabulated in Appendix 4, A 
plot of the generalized probable error values vs sharpness 
index, n is given in figure 6,12. 



CHAPTER 8 


CQECLUSIOHg 

1» A mathematical model of the coal cleaning process 
has boon derived which generates the separation (Tromp) 
curve. The 4-parametcr model reproduces the distribution 
curves with a high degree of accuracy. 

2. Rosonbrock hill climb, method with uniform weights 
to the washed coal recovery f(d) data, is suitable for 
estimating the 4 parameters in our model. 

3. The model predicts self -similar distribution 
curves in d/dg in agreement with the plant data. 

4. The model reproduces the distribution curves better 
than Gottfried's mathematical representation of distribution 
curves by the Veibull function. 

5. The estimated values of the 4 parameters have been 
analysed for all coal cleaning equipments. Higher the value 
of sharpness index, n, more efficient is .the process. For 
bett-er performance , n should- be high, f^Cd) should be maximum 
f^(d) should be minimum and d^ should be high. 

6* Except for concentrating table valuable insight 
into the performance behaviour of all other float— sink equip- 
ment has been obtained, specially with respect to the 
important coal feed sise variable. 



67 


R EPEREICES 

1. Thomas A. Hendricks, The Origin of Coal in H.H. Lowry 
(Editor), 1 Chemistry of Coal Utilizations John ¥iley 
and Sons, 1945, Ch.l. 

2. G-ilbort Theissen, Composition and Origin of Mineral 
Matter in Coal, In II, I-I, Lowry (Editor) t Chemistry of 
Coa.1 Utilization 1 , John Wiley and Sons, 1945, Ch*14. 

3. Metrics A.C., Tlio Pormation, Composition and Classifi- 
cation of Coa.1. In: P.C, Pope (Editor), 1 Coal-Production, 
Distribution, Utilization 1 , Chapman and Hall Ltd., 

1945, Ch.l. 

4. II. P. Yancey and Geer, M.R., The Cleaning of Coal In: 

II. II. Lowry (Editor), 'Chemistry ofCoal Utilization 1 , 

John Uiloy and Sons, 1945, Ch.16, p.598, 

5. Zac-har, P.E. and Gilbert A.G., The Economics of Coal 
Preparation, In: Joseph W, Leonard and David R. Mitchell 
( Editors) , 1 Coal Preparation 1 , The American Institute of 
Mining, Metallurgical and Petroleum Engineers, 1968, 
Ch.5. 

6. Major E. C. Dixon, Coke and by-products Manufacture 
Charles Griffin and Company Ltd. , 1939, Ch,2. 

7. Wilkins E.T., The Preparation of Coal for the Market, 

In: P.-C.Popo (Editor), 1 Coal-Production, Distribution, 
Utilization 1 , Chapman and Hall Ltd., 1945, Ch. 4* 



68 


8 * 


9. 


Yancey, H.P. and Goer, M.R., The Cleaning of Coal 

In: HJI. Lowry (Editor), 'Chemistry of Coal Utilization' , 
John Wiley end Sons, 1945, Ch.16, p. 572 . 

Joseph W. Leonard and David 11. Mitchell (Editors), 

‘Coal Preparation* , Tho American Institute of Mining, 

Metallurgical and Petroleum Engineers, 1968, Ch,9 and 
Ch . 10 . 


9(a) Palowitch E.R, and Deurbrouck,A.W. , Wet Concentration 
oi Coarse Coal, Part 1- Dense Medium Separation, In: 
Bef.9 page 9-3. 

9(b) Lovell, H.L. ,Wet Concentration of Coarse Coal, Part 2— 
Hydraulic Separation, In: Ref. 9 page 9-38. 

9(c) Sokashi, II., Geer, M.R. , Yoncoy, H.C*,Wet Concentration 
of Pine Coal Part 1, Dense Medium Separation, In: 

Ref. 9 page 10-3. 

9(d) D eur obrouck , A . W. and Palowitch E.R. , Wet Concentration 
of Pine Coal Part 2, Hydraulic Concentration, In: Ref. 

9 page 10,32. 

10. Gokhalc, K.V.G.K. and Rao,T.C» Coal and Lignite In: 

Oro Deposits of India. Their Distribution and Processing, 
Thomson Press (India) Ltd., 1973, Ch.18. 

11. Geor, M.R. and Yancey, H.F. , Plant Performance and 
Porecasting Cleaning Results, In: Ref. 9 Ch. 18, 

12. Byron S. Gottfried, A Generalization of Distribution 
Data for Characterizing tho Performance of Ploat-sink 



69 


CocJ. Cleaning Devices, International Journal of 
Mineral Processing, 5 (1978) 1-20. 

13 . Byron, S. Gottfried and Jacobsen, P.S.,. Generalized 
Distribution Curve for Characterising the Performance 
of Coal Cleaning Equipment. U.S. Bureau of Mines, 

Rep. Invest. 8238, 1977.' 

14. Box, M.J., 'A Mow Method of Constrained Optimization 
and a Comparison with Other Methods’ , Computer Journal 
12 , ( 1962 ). 

15. David H. Himmolblau, Nonlinear Models, In: Process 
Analysis by Statistical Methods. John Mi ley and Sons, 
Inc., 1968, Ch.6. 

16. Rosonbrock, H.H., ’An Automatic Method for Finding 
the Greatest or Least Value of a Function’ , Computer 
Journal, Vol.3, I960, page 168-, 

17. Eucstcr, J.I, and Mize, J.H., ’Optimization Techniques 
With Fortran’ , He Graw Kill, 1973, p. 386. 



70 


APPENDIX 1 

BOX-COMPLEX T EGHITIQUB [ 14—15 ] 

Phis method is Based on the flexible geometric simplex 
method* Given ‘IT 1 independent variables method requires ’ST+l* 

I » f A 

feasible starting cases (initial guesses). All the 'E+l* 
cases are randomly generated by varying the initial variable 
values between the upper and the lower bounds and are checked 
against the constraints. In this work for 4-pa-rameter models, 
5 starting cases randomly generated from the initial guess 
point* 

Once the starting cases are provided, the search phase 
begins by evaluating the objective function, In our case at 
each of 5 points, and the worst case (for which the error is 
maximum) is rejected* The centroid of the remaining cases is 
found oiit and the rejected point is reflected on the opposite 
side of the centroid by a factor of GKR, reflection coefficient 

X i,j(now) “ ® ILl ^ A i,c “ 2 i,j (rejected)/ ^i,c 
where X- . is the new point of the ith parameter 

1 f J 

i = 1, 2, .... ,21 

The value of GMB. is arbitrarily chosen as 1.1, If the 
new case violates the constraints, it is modified by either 
GI-1E, expansion coefficient or (31C, contraction coefficient, 
so that the new case is moved towards the centroid. This 



71 


procedure of moving towards the centroid is continued 
until the constraints; are no more violated. After every 
iteration the objective function is checked for convergence 
until it is achieved. Figure 7 is the flow chart for the 
program used. 
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APPENDIX 2 


PQSMBROGX HUP CL 11-33 . I- IETHOD [ 16-17] 

Unlike the Box-Complex technique this method requires 
only one feasible starting case (Initial guesses). In this 
method the initial guesses is chosen and the appropriate 
constraints are posed on all the variables. Once the feasible 
starting case is given, initial step sizes, E^(i=l, 2, . . . ,1T) 
are picked up and the objective function is evaluated. Then 
the first variable X^_ is stepped a distance 'E^, 1 parallel to 
the axis and the objective function is evaluated. The current 
best objective function is defined by a value of E 0 where 
the constraints are satisfied and , H 1 the current best 
objective function value for a point where the constraints 
arc satisfied in addition to the boundary zones are not 
violated. , H t and 1 E Q ’ arc initially set equal to the objective 
function value at the starting point. If in the current point 
objective function evaluation, II is worse than E q , or if 
the constraints arc violated, then the trial is a failure and 
the value of ’ E^, is decreased by a factor of 0w5 and the 
direction of movement is reversed. Also if there is function 
improvement then E^' is increased by three times its original 
value and the search is carried out in the next direction. 

On the otherhand, if the current point lies within a. 
boundary zone, the objective function is modified as follows: 
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Por lower zone, where Gj- is lower constraint and 
Hg is upper constraint. 


ana 


P w 


~~ ^TT ~ (Bjr-Gjr) 40 ^ ) 


E E/ 


(%-Gfc) 10' 


-4 


Por upper zone, 

At the inner edge of the boundary zone, ? w =0, that 
is the function is unaltered (P( n0TT ) = Then the 

function value is replaced by the best current function value 
in the feasible region and not in the boundary zone. Por the 
function which improves as the constraints is approached, the 
modified function has an optimum in the boundary zone. 

After search in all the coordinate directions, it is 
checked for convergence. If convergence is not achieved then 
it is chocked for one success and one failure in each 
direction and if it is not so then the whole procedure is 
repeated. On the other hand, if one success and one failure 
is achieved in each direction, eh on the axes are rotated 
according to the equation, 
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Also i = variable index (i = 1, 2, . . . . ,J) 
j = direction index (j = 1,2,..., IT) 

E = stage index 

d^ = sum of distance moved in the 1 i 5 direction since 
last rotation of axes. 

II. . = direction vector component normalized. 

1 , 3 

After rotating the axes the search is ma,de in each of 
the x direction using the new coordinate axis. 


-MO 

i(new) 


= S (K) 
X i(old) 


+ E. 


(K) 


I.T. .« 
1. J 


The above procedure is repeated until convergence is 
achieved. Figure 8 is the flow chart for the program used. 
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EFFECT OF SHARPNESS INDEX ON DISTRIBUTION CURVE » » FCD)*1 ,/(+D*»N) * 
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